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1 .  Introduction 


During  the  reporting  period  studies  have  been  conducted  in  two 
major  areas  in  fundamental  combustion,  namely  soot  formation  and 
the  structure  of  laminar  flames.  Concerning  soot  formation,  the 
emphasis  has  been  on  further  resolving  the  relative  importance  of 
temperature  versus  dilution  in  soot  formation.  As  a  companion 
project,  we  have  also  identified  the  dominant  physical  and 
chemical  influences  when  a  certain  additive,  say  carbon  dioxide, 
is  added  to  either  the  fuel  or  the  oxidizer  side  of  the  flame.  On 
flame  structure  studies,  investigations  have  been  conducted  on 
premixed  flames,  on  diffusion  flames,  and  on  the  concept  of 
flammability  limits.  For  premixed  flames,  emphasis  has  been  on  the 
determination  of  the  laminar  flame  speeds  and  the  validation  or 
identification  of  the  associated  chemical  kinetics  for  the  mixture 
of  interest.  Studies  have  been  conducted  on  all  of  the  C2  fuels  as 


well  as  hydrogen.  For  diffusion  flames,  we  have  examined  the 
accuracies  of  the  various  descriptions  of  the  flow  fields 
governing  counterflow  diffusion  flames.  The  study  further 
indicates  the  importance  of  the  strain  rate  experienced  by  the 
reaction  zone  of  the  flame  in  causing  extinction.  Concerning 
flammability  limits,  a  kinetic  criterion  has  been  proposed  from 
which  the  rich  and  lean  flammability  limits  of  a  fuel-oxidizer 
system  can  be  predicted  from  first  principles.  Agreement  between 
predictions  and  experimental  data  has  been  extremely  close, 
indicating  the  possibility  that  the  flammability  limits  are 
kinetically  controlled,  and  that  these  limits  can  be  predicted.  If 
what  we  have  found  is  indeed  correct,  after  further  testings,  it 
will  be  a  significant  milestone  in  fundamental  as  well  as  in 
practical  combustion. 


The  various  project 
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2  .  Dilution  Effects  in  Soot  Formation  in  Diffusion  Flames 


A  major  concept  developed  in  recent  studies  on  soot  formation 
in  diffusion  flames  is  that  the  flame  temperature  exerts  the 
dominant  influence  on  soot  formation.  This  concept  has  been 
experimentally  investigated  by  lowering  the  flame  temperature  by 
diluting  the  fuel  stream  with  an  inert. 

In  a  previous  study  of  ours  we  challenged  this  concept  by 
suggesting  that  fuel  dilution,  in  addition  to  the  lowering  of  the 
flame  temperature,  could  also  be  responsible  in  reducing  the 
amount  of  soot  formation,  and  that  previous  studies  were  not 
sufficiently  well  controlled  because  the  methodology  involved  the 
simultaneous  variation  in  the  fuel  concentration  as  well  as  the 
flame  temperature. 

Subsequently  we  proposed  an  alternate  methodology  through  which 
temperature  and  concentration  variations  can  be  separately 
imposed.  The  methodology  involves  controlling  the  flame 
temperature  by  substituting  a  given  amount  of  the  nitrogen  in  the 
oxidizer  by  an  equal  amount  of  argon,  which  has  a  different 
specific  heat.  Thus  the  concentrations  of  the  fuel  and  oxidizer 
streams  can  be  held  fixed  while  the  flame  temperature  is  changed. 
Consequently  we  demonstrated  that  not  only  temperature  and 
dilution  are  both  important  in  soot  formation,  in  many  instances 
dilution  can  be  even  more  important.  This  methodology  has  now  been 
termed  as  the  heat  capacity  technique. 

Our  previous  experiments  were  conducted  by  using  the 
counterflow  flame  while  most  of  the  literature  experiments  were 
performed  with  the  co-flow  flame.  Thus  reservations  have  been 
raised  concerning  the  applicability  of  our  results  to  those  of  the 
co-flow  flames.  Thus  in  this  study  we  have  applied  the  heat 
capacity  technique  to  identify  the  separate  influences  of 
temperature  and  dilution.  Figure  1  shows  the  integrated  soot 
volume  fraction  as  a  function  of  the  axial  position.  Curve  A  has 
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100%  C2H4,  Curve  B  50%C2H4/50%N2r  and  Curve  is  the  same  as  Curve  B 
except  the  flame  temperature  has  been  adjusted  to  be  equal  to  that 
of  Curve  A.  Thus  a  comparison  between  Curves  B  and  C  shows  the 
effects  due  to  temperature  variation  because  they  have  the  same 
concentration  but  different  temperatures,  while  a  comparison 
between  between  Curves  C  an  A  shows  the  effects  due  to 
concentration  variation  because  they  have  the  same  temperature  but 
different  concentrations.  The  results  are  qualitatively  similar  to 
what  was  observed  with  the  counterflow  flame,  and  dramatically 
demonstrate  the  fact  that  dilution  can  exert  a  much  stronger 
effect  than  temperature  in  soot  formation. 

This  work  is  reported  in  Publication  No.  2. 

We  mention  in  passing  that  recent  studies  by  Naegeli  at  the 
Southwest  Research  Institute  and  Gulder  in  Canada  on  premixed 
flames  have  also  substantiated  our  contention  that  dilution  is  an 
important  factor  in  soot  formation  in  flames. 


3  .  Effects  of  Additives  in  Soot  Formation 

The  use  of  gaseous  additives  has  been  frequently  suggested  as  a 
viable  strategy  to  reduce  soot  formation.  The  cause  for  such  a 
reduction,  however,  has  not  been  adequately  and  systematically 
studied.  Take  CO2  as  an  example  because  of  its  potential  for 

implementation  through  exhaust  gas  recirculation.  When  CO2  is 
added  to  a  fuel  or  oxidizer  stream,  both  the  flame  temperature  and 
the  reactant  concentration  are  reduced,  as  just  mentioned.  CO2, 
however,  can  also  participate  chemically  in  the  soot  formation 
process.  Thus  a  viable  study  of  additive  effects  must  necessarily 
identify  the  separate  effects  of  temperature  reduction,  reactant 
dilution,  and  direct  chemical  participation. 


Figure  2  shows  typical  results  obtained  for  CO2  addition  to  the 


fuel  stream  in  C2H4  flames.  Here  we  have  used  a  counterflow  flame 
and  have  plotted  the  soot  particle  inception  limit  versus  the 
additive  mole  fraction.  The  results  were  obtained  in  the  following 
manner.  First,  the  soot  inception  limits  for  various  amounts  of 
fuel  dilution  were  obtained.  The  experiments  were  repeated,  for  a 
given  amount  of  fuel  dilution,  by  replacing  a  portion  of  the 
nitrogen  in  the  oxidizer  with  an  equal  molar  amount  of  argon  so 
that  the  flame  temperature  was  that  of  the  pure  fuel.  The 
difference  then  indicates  the  thermal  effect  of  the  additive. 

We  next  replaced  CO2  by  equal  molar  amounts  of  N2,  ^nd  readjust 
the  flame  temperature  back  to  that  of  the  pure  fuel.  The 
additional  difference  observed  then  gives  the  chemical  effect. 
Finally,  the  difference  between  this  result  and  that  obtained 
without  any  dilution  yields  the  influence  of  dilution.  We  have 
therefore  succeeded  in  isolating  the  individual  effects  of 
dilution,  temperature,  and  direct  chemical  participation.  It  is 
clear  that  each  of  these  effects  accounts  for  a  measurable 
reduction  in  soot  formation  and  therefore  cannot  be  neglected. 

This  work  is  reported  in  Publication  No.  3. 


4  .  Flame  Speeds  and  Oxidation  Kinetics  of  C2~Hydrocarbons 

In  our  previous  work  for  OSR  we  have  developed  an  experimental 
method  through  which  laminar  flame  speeds  of  gaseous  combustibles 
can  be  accurately  determined.  The  method  involves  determining  the 
flame  speeds  of  strained  flames  situated  in  counterflows  whose 
strain  rates  can  be  measured.  By  plotting  the  strained  flame  speed 
versus  the  strain  rate,  the  unstrained  laminar  flame  speed  can  be 
determined  through  extrapolation  to  zero  stretch.  This  is  a 
significant  advancement  because  practically  all  previous  methods 
used  for  the  determination  of  laminar  flame  speeds  have  flame 
stretch  effects. 
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Accurate  values  of  laminar  flame  speeds  are  needed  in  many 
aspects  of  fundamental  and  practical  research.  One  recent  use  is 
the  partial  validation  and/or  determination  of  the  chemical 
kinetic  mechanisms  for  the  oxidation  of  given  fuel  systems.  It  is 
obvious  that  such  an  effort  is  meaningful  only  if  accurate 
experimental  values  of  the  laminar  flame  speeds  are  available. 

We  have  previously  determined  the  laminar  flame  speeds  and  the 
kinetic  scheme  of  methane/oxygen  mixtures  over  extensive  ranges  of 
concentration,  pressure,  and  flame  temperature.  In  the  present 
study  we  have  extended  the  investigation  to  cover  all  the  C2 

hydrocarbon  species  of  ethane,  ethylene,  and  acetylene,  as  well  as 
propane.  Such  an  extension  is  crucial  because  methane  only  has  C-H 
bonds  while  the  breaking  of  the  C2  bonds  is  important  in  the 
kinetics  of  all  of  the  higher  hydrocarbons.  Through  this  study  we 
have  found  that  all  of  the  existing  C2  kinetic  schemes  are 

inadequate  to  describe  our  experimental  result.  Consequently  we 
proposed  a  new  scheme  which  was  able  to  predict  practically  all  of 
our  experimental  data.  Figure  3  shows  the  comparison  between 
experimental  and  calculated  laminar  flame  speeds  for  ethane, 
ethylene,  and  acetylene  over  extensive  ranges  of  equivalence 
ratios  and  under  1  atmosphere  pressure.  The  agreement  can  be 
considered  to  be  quite  satisfactory. 

This  work  is  reported  in  Publication  No.  4. 


5  .  Flame  Speeds  and  Oxidation  Kinetics  of  Hydrogen 

Recent  interests  in  supersonic  combustion  have  identified  the 
importance  of  flame  stabilization  and  hydrogen  oxidation  kinetics. 
Further  discussions  with  chemical  kineticists  have  convince  us 
that  the  hydrogen  oxidation  kinetics  is  far  from  being 
established,  especially  for  the  very  weakly-burning  situations 
with  low  flame  temperatures.  This  is  certainly  disturbing  because 
hydrogen,  being  the  simplest  fuel  molecule,  forms  the  building 
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block  for  the  oxidation  of  all  hydrocarbons.  Furthermore,  the 
weakly-burning  situations  are  relevant  to  ignition,  extinction, 
flame  stabilization,  engine  knock,  etc. 

A  detailed  literature  survey  then  shows  that  there  has  not  been 
many  attempts  at  determining  the  laminar  flame  speeds  of  ultra 
lean  hydrogen/air  mixtures.  The  difficulty  here  is  that  lean 
hydrogen/air  flames  are  prone  to  dif fusional-thermal  cellular 
instability.  Consequently  smooth  flames  have  not  been  observed. 
Further  recognizing  the  highly  diffusive  nature  of  hydrogen  and 
its  strong  influence  by  stretch,  the  usefulness  of  the  previous 
data  on  lean  hydrogen-air  flame  speed  is  in  question. 

In  this  study  we  have  therefore  accurately  determined  the 
laminar  flame  speeds  of  hydrogen/air  mixtures  over  extensive 
variations  of  the  equivalence  ratio  and  pressure.  A  comparison  of 
these  data  with  calculated  values  using  different  literature 
kinetic  schemes  reveal  the  inadequacy  of  existing  mechanisms  under 
low  flame  temperatures.  Figure  4  shows  that  the  experimental 
values  are  consistently  higher  than  the  calculated  values. 
Systematic  sensitivity  analysis  and  adjustment  of  kinetic 
parameters  fail  to  improve  the  discrepancy.  It  is  therefore 
believe  that  there  is  a  missing  channel  of  chain  branching  which 
has  not  been  identified. 

This  work  is  reported  in  Publication  No.  5. 


6.  Flammability  Limits 

An  unsolved  problem  in  fundamental  combustion  research  is  the 
meaning  and  possible  existence  of  flammability  limits. 
Empirically,  it  has  been  observed  that  a  particular  fuel/oxidizer 
system  becomes  absolutely  noncombustible  when  its  (fuel) 
concentration  is  beyond  either  a  lean  or  a  rich  limit.  For 
example,  for  the  methane/air  system,  the  lean  and  rich  limits  are 
respectively  4.8  and  15%  of  methane.  The  confusion  governing 
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flammability  limits  is  that;  (a)  It  is  not  sure  whether 
flammability  limit  indeed  exists;  (b)  What  is  the  mechanism 
governing  its  onset  if  it  exists;  and  (c)  How  can  it  be  predicted 
from  first  principles? 

In  the  present  investigation  we  have  first  conceptually  defined 
the  phenomenon  of  flammability  limit  as  the  mixture  concentration 
beyond  which  steady  propagation  of  an  adiabatic,  planar,  flame 
becomes  impossible.  Clearly,  if  a  flame  cannot  propagate  in  such 
an  idealized  mode,  it  will  not  be  able  to  propagate  in  other 
nonadiabatic  situations.  Next,  in  order  to  quantitatively 
determine  a  flammability  limit  from  first  principles,  we  identify 
a  limit  as  the  state  at  which  the  termination  reaction  of  a 
crucial  chain  branching-termination  reaction  pair  becomes  more 
sensitive  to  changes  in  the  mixture  concentration  than  the 
corresponding  branching  reaction.  Thus  as  long  as  the  kinetic 
mechanism  of  the  fuel/oxidizer  system  is  known,  the  state  of 
flammability  can  be  predicted.  Finally,  we  have  also 
experimentally  determined  the  flammability  limits  of  a  great 
variety  of  mixtures  by  first  determining  the  extinction  limits  of 
an  adiabatic  stretched  flame,  and  then  extrapolating  the  results 
to  zero  stretch. 

Table  1  shows  the  presently  determined  theoretical  and 
experimental  flammability  limits  in  terms  of  the  fuel 

concentration  (Q*)  and  equivalence  ratio  (<J>*).  Figure  5  compares 
the  theoretical  and  experimental  limits.  The  closeness  of  the 
comparison  is  quite  evident.  The  study  further  identifies  many 
properties  of  the  flammability  limits. 

This  work  is  reported  in  Publication  No.  6. 


7  .  Extinction  of  Methane/Air  Diffusion  Flames  by  Stretch 

Although  extensive  amount  of  work  has  been  performed  on  the 
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structure  and  extinction  of  laminar  diffusion  flames,  there  does 
not  exist  any  quantitative  experimental  data  on  the  flame 
structure  determined  by  using  laser  techniques.  Thus  in  the 
present  investigation  we  have  first  accurately  determined  the 
velocity  structure  of  counterflow  diffusion  flames  by  using  LDV, 
especially  for  states  close  to  extinction.  These  profiles  are  then 
compared  with  numerically-calculated  ones  with  different 
descriptions  of  the  flow  field.  Figure  6  shows  that  the  comparison 
is  quite  good  if  the  flow  field  is  described  by  the  plug  flow.  It 
is  further  found  that  the  maximum  velocity  gradient,  which  occurs 
on  the  fuel  side  of  the  main  reaction  zone,  achieves  a  value  at 
extinction  that  is  relatively  insensitive  to  the  boundary 
conditions  of  the  flow. 

This  work  is  reported  in  Publication  No.  7. 
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Figure  6 


16 


TABLE  1 


(L/R)  Mixture 

P. 

atm 

Experimental 

Tad, 

K 

S° 

cm/s 

Theoretica] 

(L) 

CHVair 

3.00 

5.18  (0.52) 

1520 

— - 

5.46  (0.55) 

(L) 

CHi/air 

2.00 

5.08  (0.51) 

1500 

.... 

5.18  (0.52) 

(R) 

CHVair 

2.00 

15.5  (1.75) 

1727 

4.0 

16.1  (1.82) 

(L) 

Ciii/air 

1.00 

4.80  (0.48) 

1442 

4.3 

4.80  (0.48) 

(R) 

CHi/air 

1.00 

14.9  (1.67) 

1782 

5.3 

15.1  (1.69) 

(L) 

CHii/air 

0.75 

4.90  (0.49) 

1461 

5.0 

4.90  (0.49) 

(L) 

CRi/air 

0.50 

4.80  (0.48) 

1481 

6.0 

4.80  (0.48) 

(L) 

CH^air 

0.375 

5.08  (0.51) 

1480 

9.0 

4.51  (0.45) 

(L) 

CH4/air 

0.25 

4.99  (0.50) 

1470 

11.1  4.32  (0.43) 

(L) 

CH4/"Ar-air " 

1.0 

3.45  (0.34) 

1412 

3.9 

3.55  (0.35) 

(R) 

CH4/"Ar-air " 

1.0 

17.9  (2.07) 

1771 

7.5 

17.4  (2.00) 

(L) 

CH4/i0.5925 Ar+0.  i  975He+0V2 1 62]^ 

1.0 

3.84  (0.38) 

1525 

10.4  3.55  (0.35) 

(L) 

CHVi0.395Ar+6.395He-<^‘.2 1 62]* 

1.0 

4.23  (0.42) 

1634 

20.5  3.60  (0.36) 

(L) 

CiviO.  1975Ar+b.5925He+0.2i  62]“ 

1.0 

4.42  (0.44) 

1687 

31.5  3.55  (0.35) 

(L) 

CH4/"He-aii  " 

1.0 

4.70  (0.47) 

1765 

43.8  3.55  (0.35) 

(L) 

C^[b.i802+6.82N2]' 

1.0 

4.96  (0.58) 

1477 

4.7 

4.88  (0.57) 

(L) 

CH4/'[b.iM02+b.849N2]® 

1.0 

4.95  (0.69) 

1475 

5.0 

4.81  (0.67) 

(R) 

CH4/ib.iM02+b.849N2]“ 

1.0 

9.06  (1.32) 

1682 

4.6 

9.25  (1.35) 

(L) 

i6.94CH4-t^.66CH3Bri/air® 

1.0 

5.42  (0.57) 

1597 

— - 

5.15  (0.54) 

(L) 

H2/air 

1.5 

6.52  (0.17) 

828 

.... 

—  ( — ) 

(L) 

ii^air 

1.0 

6.30  (0.16) 

811 

.... 

....  (....) 

(L) 

H2/air 

0.5 

5.74  (0.14) 

767 

.... 

—  ( — ) 

(R) 

Hi/0^2  (^‘=1.6? 

1.0 

15.2  (--) 

1067 

7.0 

13.1  (--) 

(R) 

Hi/02/N2  (^'=1.6? 

0.75 

11.9  (-.-) 

903 

1.0 

11.1  (--) 

(R) 

Hi/b2/N2  (<(>=1.6)^ 

0.5 

8.49  (--) 

736 

0.5 

7.43  (--) 

(R) 

Hi/02/N2 

1.0 

19.8  (--) 

1101 

12.0  16.9  (— -) 

(L) 

(a25CH4'-1^.75H2i/air* 

1.0 

5.66  (0.25) 

998 

.... 

....  (....) 

(L) 

(6.56CH4+b.5bH2i/air‘ 

1.0 

5.10  (0.32) 

1136 

.... 

....  (....) 

(L) 

i6.75CH4-^.25H2i/air‘ 

1.0 

5.03  (0.41) 

1313 

.... 

5.27  (0.43) 

(L) 

C2H6/air 

1.0 

2.91  (0.50) 

1506 

7.2 

2.80  (0.48) 

(R) 

C2H6/[b.ii^2^.’84N2]‘ 

1.0 

7.59  (1.80) 

1512 

3.0 

7.20  (1.70) 

(L) 

QlVair 

1.0 

2.92  (0.43) 

1516 

5.8 

2.79  (6.41) 

(R) 

C2H4/[b.  i'202+6.88N2]^ 

1.0 

6.44  (1.72) 

1484 

3.2 

6.09  (1.62) 

(L) 

C2H2/air 

1.0 

2.86  (0.35) 

1400 

9.0 

2.78  (0.34) 

(L) 

CsHg/air 

1.0 

2.06  (0.50) 

1511 

7.6 

1.94  (0.47) 

(L) 

[0.92CO+0.08H2]/air 

1.0 

11.2  (0.30) 

1277 

1.8 

11.9  (6.32) 

(L) 

ib.92CO-K).08H2]/[0. 1 1 3O2+0.887N2i* 

1.0 

13.0  (0.66) 

1408 

4.1 

13.0  (0.66) 

(R) 

(6.92Cb-^.68H2j4b.  1  i  302H^.887N2i‘ 

1.0 

52.5  (4.90) 

1245 

8.5 

54.0  (5.20) 

•  Species  coefficients  represent  mole  fractions  in  fuel  and/or  oxidizer  mixtures. 
^  Flammability  limit  obtained  by  varying  the  nitrogen  dilution  for  fixed  <|>. 
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COMBUSTION  AND  FLAME  80;  7- 16  ( 1990)  7 

Chain  Mechanisms  in  the  Overall  Reaction  Orders  in  Laminar 

Flame  Propagation 

F.  N.  EGOLFOPOULOS  and  C.  K.  LAW 


Department  of  Mechanical  and  Aerospace  Engineering,  Princeton  University,  Princeton,  NJ  0S544 


The  laminar  flame  speeds  of  methane-oxygen-nitrogen  mixtures  as  functions  of  the  flame  temperature  and 
system  pressure  p  have  been  experimentally  determined  by  using  the  counierflow,  twin-flame  technique.  These 
data  are  then  compared  with  numerically-calculated  values  obtained  by  using  an  independently  validated  kinetic 
scheme.  Results  show  that  the  overall  reaction  order  n  and  activation  energy  are  far  from  being  constants, 
that  n  decreases  with  increasing  p  and  decreasing  whereas  E,  increases  with  increasing  p.  and  that  n  can 
actually  assume  negative  values;  the  last  result  implies  that  the  mass  burning  rate  of  some  weakly  burning  flames 
may  decrease  with  increasing  pressure.  By  further  identifying  the  crucial  reaction  steps  through  sensitivity  analysis, 
the  present  results  are  interpreted  on  the  basis  of  the  influence  of  chain  branching-termination  mechanisms  on  the 
overall  reaction  rate. 


INTRODUCTION 

Chain  mechanisms,  especially  those  involving  the 
competition  between  branching  and  termination 
reactions,  are  of  intrinsic  importance  to  chemi¬ 
cally  reacting  systems  because  they  are  potentially 
capable  of  qualitatively  modifying  the  combustion 
behavior  from  those  anticipated  on  the  basis  of, 
say,  one-step  overall  reactions.  Examples  pertain¬ 
ing  to  purely  chemical  systems  include  the  explo¬ 
sion  limits  of  hydrogen-oxygen  mixtures  and  the 
cool  flame  phenomena  of  hydrocarbon-air  mix¬ 
tures.  Comparatively  fewer  studies,  however,  have 
been  conducted  on  the  influence  of  branching- 
termination  chain  mechanisms  on  the  global  be¬ 
havior  of  /lames,  for  which  molecular  transport 
as  well  as  chemical  kinetic  processes  are  inher¬ 
ently  coupled.  Recently  the  interest  and  activities 
in  this  class  of  phenomena  have  increased  con¬ 
siderably  because  of  the  increased  capability  to 
conduct  detailed  numerical  modeling  of  combus¬ 
tion  phenomena.  The  particular  flame  configura¬ 
tion  that  has  been  frequently  adopted  for  such  stud¬ 
ies  is  the  steady  propagation  of  the  adiabatic  one¬ 
dimensional  planar  flame  in  the  doubly-infmite  iu- 
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main.  The  flow  field  in  this  configuration  is  well 
defined  and  is  the  simplest  possible,  thus  allow¬ 
ing  the  unambiguous  exploration  and  identifica¬ 
tion  of  the  influence  of  chemical  kinetics  on  the 
flame  behavior.  Consequently,  through  numerical 
simulation  of  such  flames  Westbrook  and  Dryer 
[1,  2],  Wamatz  [3],  Dixon-Lewis  [4],  and  Coffee 
and  co-workers  [5,  6]  have  amply  demonstrated 
the  importance  of  chain  mechanisms  in  the  flame 
structure  and  propagation.  It  is  further  emphasized 
[1,2]  that,  because  of  chain  termination,  the  cal¬ 
culated  laminar  flame  speed  decreases  with  pres¬ 
sure  at  a  faster  rate  when  the  pressure  exceeds 
a  certain  range.  Because  this  behavior  is  central 
to  the  present  study,  we  briefly  discuss  the  cause 
responsible  for  it. 

First,  a  simple  heat  and  mass  balance  shows 
that  the  burning  rate  (per  unit  area)  of  the  lam¬ 
inar  flame,  /n°,  is  related  to  the  overall  reaction 
rate  w  through  ~  Since  for  a  given  mix¬ 
ture  concentration  and  upstream  conditions  we  can 
write 

w  ~ p"  exp( -Ea /R°T ti),  (I ) 
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where  p  is  the  pressure,  the  adiabatic  flame 
temperature,  and  n  and  Eg ,  respectively,  the  over¬ 
all  reaction  order  and  activation  energy  pertaining 
to  that  particular  mixture  and  burning  configura¬ 
tion,  we  have 

-p"- exp(-£„/2/?‘>r.d).  (2) 

In  writing  Eq.  2  it  is  recognized  that  the  re¬ 
maining  influences  of  pressure  and  flame  temper¬ 
ature,  through  transport  property  variations,  are 
secondary. 

The  parameters  n  and  Eg  are  to  be  empirically 
determined.  Thus  with  a  given  set  of  either  ex¬ 
perimental  or  numerical  data  m^{p,  T^a),  these 
parameters  can  be  defined  as 


The  mass  burning  rate  is  the  eigenvalue 
of  laminar  flame  propagation  and  therefore  is 
the  proper  parameter  to  indicate  the  burning  in¬ 
tensity  of  the  flame.  However,  traditionally  the 
laminar  flame  speed  has  been  used.  Since 
Eq.  2  can  be  expressed  in  terms  of 
as 


Si  ^py  cxp(-Eg/2R°Tg^),  (5) 

where  pg  is  the  density  of  the  unburned  mixture 
and 


Equation  6  shows  that  the  dependence  of  Si  on 
p  includes  the  effects  of  both  chemical  reaction 
(n/2)  and  density  variation  (-1).  Thus  when  Sj 
is  used  for  data  interpretation,  the  density  effect 
has  to  be  first  subtracted  out  before  inferences 
can  be  made  regarding  the  effect  of  pressure  on 
the  burning  intensity. 

We  next  consider  the  influence  of  chain  mecha¬ 


nisms  on  the  laminar  burning  rate  via  the  example 
H  -I-  O2  —  OH  -I-  O,  (B) 


H -1-02 -l-M HO2 -t-M,  (T) 

which  are  two  of  the  crucial  chain  branching  and 
(quasi-)  termination  reaction  steps  in  the  oxida¬ 
tion  kinetics  of  hydrogen  and  hydrocarbons.  The 
branching  step  (B)  is  a  temperature-sensitive,  two- 
body  reaction,  while  the  termination  step  (T) 
is  a  temperature-insensitive,  three-body  reaction. 
Thus  the  branching  and  termination  reactions  can 
be  enhanced  relative  to  each  other  by  increasing 
the  flame  temperature  (Tad)  and  system  pressure 
(p),  respectively.  This  influence  should  then  be 
reflected  in  the  overall  reaction  rate  in  that  it 
is  increased  with  increasing  but  is  decreased 
with  increasing  p.  This  is  the  reasoning  of  West¬ 
brook  and  Dryer  [1,2]  for  their  numerical  re¬ 
sult  of  decreasing  Su°  with  increasing  p.  It  also 
supports  previous  experimental  observations:  for 
example,  Lewis  [7]  reported  a  maximum  nega¬ 
tive  value  of  7  of  about  -0.4  for  weakly  burning 
flames  whereas  Andrews  and  Bradley  [8]  showed 
that  for  methane-air  mixtures  7  =  -0.5  for  p  >  5 
atm. 

Further  insight  into  the  influence  of  chain  mech¬ 
anisms  can  be  gained  from  these  previous  re¬ 
sults.  First,  since  Sl  is  not  the  proper  eigen¬ 
value  of  laminar  flame  propagation  when  pres¬ 
sure  variation  is  involved,  a  decreasing  trend  of  Sl 
with  pressure,  OSl/dp  <  0,  does  not  provide  di¬ 
rect  indication  of  a  negative  kinetic  influence,  due 
to  chain  termination,  on  the  burning  rate  of  the 
flame.  Indeed,  by  u..lng  Eq.  6,  the  previously  re¬ 
ported  7  all  correspond  to  positive  values  of  n, 
with  2  >  n  >  1.  On  the  other  hand,  the  facts  that 
these  n  values  are  smaller  than  2,  t'at  n  also  de¬ 
creases  with  pressure,  and  that  removal  of  the 
termination  reaction  (T)  changes  the  value  of  n, 
as  pointed  out  by  Westbrook  and  Dryer  [1,  2], 
do  indicate  the  progressive  importance  of  termi¬ 
nation  reactions  with  pressure.  That  is,  since  a 
strongly  burning  flame,  dominated  by  two-body 
branching  and  carrying  reactions,  should  have  an 
n  close  to  2,  results  of  n  substantially  smaller 
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than  2  indicate  the  possible,  simultaneous  impor¬ 
tance  of  the  three-body  termination  reactions  that 
generally  tend  to  reduce  n,  and  of  the  first-order 
unimolecular  decomposition  reactions  as  the  high- 
pressure  limit  of  the  Lindemann  mechanism.  The 
relative  importance  of  the  unimolecular  reactions 
as  compared  to  the  termination  reactions  can  be 
further  eliminated  because,  if  they  were  impor¬ 
tant.  n  would  increase  with  pressure,  indicating 
the  progressive  dominance  of  the  two-body  prop¬ 
agation  reactions.  Thus  it  can  be  concluded  that 
the  previous  results  of  2  >  n  >  1  and  n  decreas¬ 
ing  with  pressure  do  provide  strong  evidence  of 
the  importance  of  the  three-body  termination  re¬ 
actions  in  the  flame  burning  intensity. 

The  present  study  aims  to  extend  these  previ¬ 
ous,  worthwhile  contributions  along  the  following 
directions.  First,  in  view  of  the  recent  recogni¬ 
tion  that  previous  experimental  determination  of 
the  laminar  flame  speeds  could  have  been  compli¬ 
cated  by  stretch  effects  [9],  we  shall  perform  ac¬ 
curate  determinations  of  8°  of  methane-02-N2 
mixtures  as  functions  of  pressure  by  using  the  re¬ 
cently  developed  counterflow  technique  [10-14]  in 
which  stretch  effects  are  systematically  subtracted 
out.  The  increased  accuracy  is  especially  impor¬ 
tant  for  the  slow-burning  flames  of  interest  in  this 
study,  as  is  further  discussed  in  the  next  section. 

We  shall  also  determine  S2  as  function  of  the 
flame  temperature,  which  yields  values  of  the 
overall  activation  energy.  Very  little  accurate  ex¬ 
perimental  data  exist  for  this  parameter. 

In  addition  to  the  experimental  investigations, 
we  shall  perform  numerical  simulations  of  these 
flames  by  using  a  kinetic  scheme  independently 
validated  for  flame  speed  calculations,  and  com¬ 
pare  the  calculated  values  with  our  experimental 
data.  Very  close  agreement  is  obtained,  thereby 
providing  confidence  in  the  accuracy  of  the  present 
study  as  well  as  the  meaningfulness  of  the  insights 
gained  herein. 

An  auxiliary  objective  of  the  present  study  is  to 
explore  the  possible  existence  of  reaction  regimes 
in  which  the  termination  reactions  are  so  strong 
that  the  overall  reaction  orders  may  assume  val¬ 
ues  of  /J  <  1  or  even  n  <  0.  The  need  for  such  an 
exploration  is  that  it  is  not  clear  a  priori  whether 
suvh  weakly  burning  flames  can  indeed  exist.  If 


they  do  not  exist,  then  the  range  of  possible  val¬ 
ues  of  n  is  significantly  narrowed.  This  can  lead  to 
considerable  simplifications  in  the  extent  of  inves¬ 
tigation  in  future  parametric  and  modeling  studies. 
On  the  other  hand,  if  they  do  exist,  especially  the 
n  <0  situations  for  which  the  mass  burning  rate 
decreases  with  increasing  pressure,  then  their  sig¬ 
nificance  should  be  recognized  and  accounted  for 
in  future  efforts. 

The  experimental  and  numerical  aspects  of  the 
present  investigation  are  respectively  presented  in 
the  next  two  sections. 

EXPERIMENTAL  METHODOLOGY  AND 
RESULTS 

The  experiment  adopts  the  counterflow  technique 
of  Law  and  co-workers  [10-14]  in  determining  the 
flame  speed.  It  basically  involves  the  establish¬ 
ment  of  two  symmetrical,  planar,  nearly  adiabatic 
flames  in  a  nozzle-generated  counterflow  config¬ 
uration,  and  determining  the  axial  velocity  profile 
along  the  centerline  of  the  flow  by  using  LDV. 
The  minimum  point  of  the  velocity  profile  can 
be  identified  as  a  reference  upstream  flame  speed 
Su,  with  the  velocity  gradient  ahead  of  it  as  the 
imposed  stretch  rate  K.  Thus  by  plotting  Sy  ver¬ 
sus  K,  we  can  obtain  the  stretch-free  flame  speed 
S2  through  linear  extrapolation  to  A  =0.  The 
entire  burner  is  housed  in  a  large  strainless-steel 
chamber  with  continuous  air  ventilation,  thereby 
allowing  for  pressure  adjustments.  Details  of  the 
experimental  setup  can  be  found  in  Refs.  12  and 
13. 

Accuracy  of  the  present  methodology  is  mainly 
limited  by  the  uncertainty  associated  with  veloc¬ 
ity  measurements,  and  is  typically  between  0.5 
and  1.5  cm/s.  This  uncertainty  range  is  generally 
smaller  than,  or  of  the  same  magnitude  as,  the 
size  of  the  data  symbols  shown  in  the  following 
figures.  Thus  error  bars  are  not  shown. 

This  methodology  is  well-suited  for  the  present 
investigation  of  pressure  effects  on  weakly  burning 
flames.  This  is  because  the  conventional  method  in 
studying  pressure  effects  on  flames  speeds  is  that 
of  the  expanding  spherical  flame  in  an  enclosed 
bomb.  For  these  spherical  flames  it  is  well  recog¬ 
nized  that  experimental  accuracy  is  low  during  the 
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initial  period  when  effects  of  flame  curvature  and 
heat  loss  to  the  electrodes  are  significant.  Further¬ 
more,  as  the  flame  grows  larger  and  data  can  be 
normally  taken,  buoyancy  will  have  caused  severe 
distortion  and  upward  displacement  of  the  flame 
ball  configuration  for  weakly  burning  flames.  Heat 
loss  to  the  electrodes  may  still  be  important  at  this 
stage  because  these  flames  are  extremely  weak.  It 
is  therefore  clear  that  accurate  flame  speed  data 
cannot  be  readily  determined  or  unambiguously 
extracted  by  using  the  spherical  flame  methodol¬ 
ogy- 

For  the  present  counterflow  methodology,  up¬ 
stream  heat  loss  is  minimized  because  the  flow  is 
nozzle-generated  and  downstream  heat  loss  is  min¬ 
imized  due  to  symmetry.  The  buoyancy  effects  are 
automatically  included  in  that  the  LDV-velocity 
profile  is  the  local  value  experienced  by  the  spe¬ 
cific  flame  segment  (around  the  stagnation  point) 
under  measurement;  in  the  Appendix  we  provide 
further  discussion  and  experimental  data  that  sub¬ 
stantiate  this  statement.  This  methodology  is  there¬ 
fore  ideally  suited  for  the  present  experiment. 

Methane-oxygen-nitrogen  flames  will  be  stud¬ 
ied  because  this  is  the  only  hydrocarbon  fuel 
whose  oxidation  scheme  has  been  extensively  in¬ 
vestigated  and  reasonably  well  established  as  far 
as  the  flame  speed  calculation  is  concerned.  Thus 
there  is  reasonable  confidence  in  comparing  the 
experimental  data  with  the  numerical  calculations 
for  these  flames. 

Two  series  of  experiments  were  conducted.  The 
first  involved  methane-air  mixtures  with  different 
equivalence  ratios,  <i),  and  at  different  pressures. 
Thus  for  a  given  0,  differentiation  of  with  p 
according  to  Eq.  3  yields  n.  A  similar  determina¬ 
tion  of  Ea  by  differentiating  by  r„i  according 
to  Eq.  4,  however,  is  less  meaningful  because  of 
the  simultaneous  influence  of  <t)  and  r,<j  on  m°. 
This  ambiguity  is  removed  in  the  second  series 
of  experiments  in  which  the  methane-oxygen  sto¬ 
ichiometry  is  fixed,  at  <^)  =  1,  and  r,<j  is  varied  by 
changing  the  nitrogen  concentration.  Thus  by  also 
varying  the  pressure,  n  and  Ea  can  be  determined 
according  to  Eqs.  3  and  4.  Calculated  7",^  for  the 
various  mixtures  are  given  in  Tables  1  and  2. 

Figure  1  shows  the  experimentally  determined 
laminar  flame  speeds  of  stoichiometric  and  lean 
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TABLE  1 

Calculated  Adiabatic  Flame  Temperatures  (K)  of  Methane- Air 
Mixtures 


Equivalence  Ratio  (<t>) 

P(atm) 

0.6 

0.8 

10 

1.2 

14 

0.25 

1668 

1996 

2197 

2130 

1977 

0.5 

1668 

2000 

2215 

2135 

1979 

1.0 

1669 

2003 

2232 

2138 

1980 

1.5 

1669 

2005 

2241 

2139 

1980 

2.0 

1669 

2006 

2247 

2140 

1980 

2.5 

1669 

2007 

2251 

2141 

1981 

3.0 

1669 

2007 

2254 

2141 

1981 

3.5 

1669 

2008 

2257 

2142 

1981 

4.0 

1669 

2008 

2260 

2142 

1981 

4.5 

1669 

2008 

2262 

2142 

1981 

methane-air  mixtures  with  different  equivalence 
ratios  and  pressures.  The  results  qualitatively 
agree  with  the  numerical  calculations  of  West¬ 
brook  and  Dryer  [1 , 2]  in  that  S2  decreases  mono- 
tonically  with  increasing  pressure. 

In  Fig.  2a  the  S2  is  converted  to  the  laminar 
burning  rate  nf  =  which,  as  discussed 

previously,  is  a  more  direct  parameter  represent¬ 
ing  the  rate  of  the  overall  chemical  reaction.  It  is 
seen  that  now  increases  with  pressure  for  all 
values  of  <t>.  By  fitting  these  data  and  then  differ¬ 
entiating  the  fitted  results  according  to  Eq  3.  the 
overall  reaction  orders  n  have  been  computed  and 
are  shown  in  Fig.  2b.  The  fitting  was  performed 

TABLE  2 

Calculaled  Adiabatic  Flame  Temperatures  (K)  of  Stoichiome¬ 
tric  Methane-Oxygen-Nitrogen  Mixtures 


N;  (%) 


P(aim) 

79 

82 

83 

84 

85  6 

1.0 

2232 

2069 

2004 

1936 

1819 

1.5 

2241 

2074 

2007 

1938 

1820 

2.0 

2247 

2077 

2009 

1939 

1820 

2.5 

2251 

2079 

2010 

1940 

1821 

3.0 

2254 

2080 

2012 

1941 

1821 

3.5 

2257 

2082 

2012 

1942 

1821 

4.0 

2260 

2083 

2013 

1942 

1822 

4.5 

2262 

2084 

2014 

1943 

1822 

21 
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flame  speeds,  ,  of  lean  methane-air  mixtures  as  functions 
of  stoichiometry  and  pressure. 

by  using  either  graphical  interpolation  or  second- 
and  third-order  polynomial  and  logarithmic  fiinc- 
tions.  The  results  demonstrate  the  following  points 
regarding  the  relative  influence  of  termination  ver¬ 
sus  branching  reactions  in  the  chain  mechanism: 


0  12  3  4 


Pressure,  atm 

Fig.  2.  Experimentally  and  numerically  determined  (a)  mass 
burning  rates,  m”,  and  (b)  overall  reaction  orders,  n,  of  lean 
methane-air  mixtures  as  functions  of  stoichiometry  and  pres¬ 
sure. 


(1)  Termination  reactions  are  enhanced  with  in¬ 
creasing  pressure,  hence  the  decreasing  trend  of 
n  with  p  for  fixed  <t>.  (2)  Branching  reactions  are 
weakened  with  decreasing  (t),  hence  the  decreasing 
n  with  decreasing  0  for  fixed  p.  (3)  n  approaches 
2  only  for  very  low  pressures,  indicating  the  pos¬ 
sibility  that  «  w  2  could  be  the  maximum  value  of 
the  reaction  order  for  this  class  of  mixtures,  for 
which  the  two-body  branching  and  carrying  reac¬ 
tions  represent  the  highest  molecularity  of  those 
elementary  reactions  that  have  positive  influence 
on  the  overall  reaction  rate.  (4)  In  view  of  (3),  it 
is  reasonable  that  with  increasing  pressure  n  de¬ 
creases  rapidly  from  2  such  that  it  is  already  close 
to  unity  even  for  the  normally  robustly  burning 
stoichiometric  mixture  at  1  atm.  (5)  The  decrease 
in  n,  however,  slows  down  after  a  certain  pres¬ 
sure  range.  This  agrees  with  the  previous  result 
18]  of  a  constant  n(  =  l)  for  p  >  5  atm  for  the  sto¬ 
ichiometric  mixture.  The  slightly  different  value  of 
the  present  n,  when  extrapolated  to  5  atm,  is  not 
surprising  considering  the  uncertainties  associated 
with  previous  measurements.  (6)  This  set  of  data 
does  not  show  negative  values  of  n.  although  n 
does  become  quite  close  to  zero  for  the  very  lean 
situations  of  0  =  0.6. 

Figures  3  and  4  show  the  corresponding  data  for 
rich  methane-air  mixtures.  The  general  behavior 
is  similar  to  those  of  the  lean  mixtures,  except 
that  n  appears  to  be  somewhat  insensitive  to  the 
stoichiometry,  at  least  for  the  two  cases  studied. 

To  further  explore  the  possible  existence  of  neg- 


Fig.  3.  Experimentally  and  numerically  determined  laminar 
(lame  speeds,  ,  of  rich  methane-air  mixtures  as  functions  of 
stoichiometry  and  pressure. 
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Fig.  4  Experimenully  and  numerically  determined  (a)  mass 
burning  rates,  m”,  and  (b)  overall  reaction  orders,  n,  of  rich 
methane-air  mtxtures  as  functions  of  stoichiometry  and  pres¬ 
sure. 

alive  n,  Figs.  5  and  6  show  variations  of 
and  n  with  pressure  for  stoichiometric  mixtures 
with  different  amounts  of  nitrogen  dilution  in  the 
(N2  +02)  mixture.  Figure  6  demonstrates  clearly 
that  although  n  remains  positive  for  slightly  di¬ 
luted  mixtures,  it  can  indeed  assume  negative  val¬ 
ues  for  the  highly  diluted  mixtures  as  pressure  in¬ 
creases. 

The  data  of  Fig.  6a  are  replotted  in  Fig.  7  to 
show  the  dependence  of  on  T^.  It  is  seen  that 
for  a  given  pressure  an  approximate  Arrhenius  re¬ 
lation  exists  over  the  broad  adiabatic  flame  tem¬ 
perature  range  of  about  1820  to  2250  K  (Table  2). 
From  these  relations  the  overall  activation  energies 
can  be  plotted  as  a  function  of  pressure  for  the  sto¬ 
ichiometric  mixture,  as  shown  in  Fig.  8.  It  is  seen 
that  Eg  increases  from  about  51.5  kcal/mol  at  I 
atm  to  75.1  kcal/mol  at  3  atm.  Recognizing  the 
progressive  importance  of  the  termination  reac- 
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Fig.  5.  Experimentally  and  numerically  determined  lami¬ 
nar  flame  speeds,  S°,  of  diluted  stoichiometric  methane- 
oxygen-nitrogen  mixtures  as  functions  of  stoichiometry  and 
pressure. 

tions  and  thereby  the  increased  difficulty  for  reac¬ 
tions  to  proceed  as  pressure  increases,  the  increas¬ 
ing  trend  of  Eg  is  reasonable.  It  is  also  significant 
to  note  the  very  large  values  of  Eg  for  the  weakly 
burning  mixtures  under  elevated  pressures. 


burning  rates,  m°,  and  (b)  overall  reaction  orders,  n,  of  diluted 
stoichiometric  methane-oxygen-nitrogen  mixtures  as  functions 
of  inert  dilution  and  pressure. 
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Fig.  7.  Arrhenius  plot  of  the  laminar  burning  rate  of  the  di¬ 
luted,  stoichiometric  methane-oxygen-nitrogen  flames. 


NUMERICAL  SIMULATION  AND 
COMPARISON  WITH  EXPERIMENTAL 
DATA 

In  order  to  calculate  the  flame  speed  and  per¬ 
form  the  sensitivity  analysis,  we  used  the  flame 
code  of  Refs.  15  and  16  with  a  small  improve¬ 
ment  in  the  computation  of  the  mass  diffusion  co¬ 
efficient  [13],  The  kinetic  scheme  used  is  the  C|- 
mechanism  of  Kee  et  al.  [16],  with  17  species  and 
58  reactions.  This  scheme  has  been  extensively 
tested  and  validated,  over  wide  ranges  of  stoi- 
chiometr)',  pressure,  and  flame  temperature  [14], 
for  flame  speed  calculations  of  lean  to  slightly  rich 
methane-oxygen-inert  mixtures.  Selected  calcula¬ 
tions  have  also  been  performed  by  using  the  full 


C2-mechanism  of  Warnatz  [17];  the  results  show 
practically  no  influence  on  the  calculated  flame 
speeds  as  well  as  the  understanding  and  conclu¬ 
sions  drawn  herein. 

The  calculated  results  are  compared  with  al¬ 
most  all  of  the  experimental  data  of  Figs.  1-8;  the 
only  exceptions  are  those  for  which  numerical  con¬ 
vergence  becomes  excessively  difficult.  The  close 
agreement  in  the  comparison  is  apparent.  Of  par¬ 
ticular  significance  is  the  result  that  negative  over¬ 
all  reaction  orders  are  also  exhibited  by  the  numer¬ 
ical  simulation.  Because  the  numerical  code  and 
kinetic  scheme  have  been  previously  validated  for 
flame  speed  calculations,  such  close  agreements, 
over  the  extensive  range  of  parametric  variations, 
strongly  support  the  accuracy  of  the  experimental 
data  and  thereby  the  viability  of  the  present  find¬ 
ing. 

From  the  calculated  temperature  profiles  the 
laminar  flame  thickness  6®  can  also  be  determined. 
Adopting  the  definition  of  Ref.  18,  the  calculated 
6®  for  the  d>  =  1.  diluted  flames  are  shown  in  Fig. 
9.  Comparing  with  Fig.  6a,  the  results  demon¬ 
strate  the  anticipated  relation 


Specifically,  for  the  highly  diluted  situations  it  is 
seen  that  6°  first  decreases  and  then  increases  with 
p,  just  as  rrP  first  increases  and  then  decreases 
with  p. 

To  further  identify  the  specific  elementary  re- 


Fig.  8.  Experimentally  determined  overall  activation  energies 
of  diluted,  stoichiometric  methane-oxygen-nitrogen  flames. 


Pressure,  atm 


Fig  9.  Numerically  determined  laminar  flame  thicknesses.  i“. 
of  diluted,  stoichiometric  methane-oxygen-nitrogen  flames 
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Fig.  10.  Normalized  sensitivities  on  mass  burning  rate  for  sto¬ 
ichiometric  methane-air  flames  at  1  and  4  atm. 

actions  that  are  responsible  for  the  observed  be¬ 
havior,  Figs.  10  and  1 1  show  the  calculated  sensi¬ 
tivities  [15,  16]  of  the  flame  burning  rates  on  the 
various  elementary  reactions  for  the  CH4-O2-N2. 
<t>  =  1  mixtures,  witli  N2  =  79  and  84%  and 
p  =  1  and  4  atm.  The  results  show  that  the 
dominant  branching  reaction  is  (B)  as  anticipated. 
There  are  two  termination  reactions  of  approxi¬ 
mately  equal  importance,  namely  reaction  (7)  and 

CHj  -h  H  M  —  CH4  -i-  M,  (T') 

which  eliminates  both  the  CHj  and  the  H  radicals. 


H+0+M=OH+M  1  1  atm 

HC(W32=H02.C0  ^  ^ 

CH4■^OH=CH3-^H20 

n 

CH4-rH=CH3-rH2 

■ 

HC0+H=C0+H2 

CH3+H+M=CH4-rM 

H-K32-rM=H02-t-M 

,  0H+H2=H20-rH 

CH4/C^ /Nj 

CH3-rOH=CH2+H20 

N2  =  84% 

0  CH3-K)=CH20+H 

0=1 

^  HCO+M=H-rCCkM 

Ba  CCHOH=C02+H 

H-*O2=0H+O 
*  •  ' 

-0.4  -0.2  0.0  0.2  0.4  0.6  0.8 

Normalized  Sensitivity  on  Mass  Burning  Rate 

Fig.  1 1 .  Normalized  sensitivities  on  mass  burning  rate  for  di¬ 
luted  stoichiometric  methane-oxygen-nitrogen  flames  at  1  and 
4  atm 
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HCO+H=CO+H2 

CH3-rH-fM=CH4-rM 
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Fig.  12.  Normalized  sensitivities  on  mass  burning  rate  for 
atmospheric-pressure  methane-air  flames  with  =  0.6  and 
1.0. 

Figure  12  shows  the  sensitivities  for  the  1- 
atm  methane-air  flames  with  <i>  =  0.6  and  1.0. 
It  is  seen  that  although  (T)  and  (T')  are  both 
important  at  0  =  1.0,  only  (T)  remains  to  be 
the  dominant  termination  reaction  at  0  =  0.6. 
This  explains  the  result  that  although  n  remains 
positive  for  the  0  =  0.6,  methane-air  flame,  it 
assumes  negative  values  for  the  0  =  1,  84/^- 
diluted  methane-Oi-N:  flames  even  though  the 
latter  have  higher  flame  temperatures.  That  is.  in 
the  former  case  there  is  only  one  dominant  termi¬ 
nation  reaction  to  slow  down  the  overall  reaction 
rate,  whereas  in  the  latter  case  there  are  two  of 
them. 

Finally,  in  Figs.  13  and  14  we  have  plotted  the 


Fig.  13.  Numerically  calculated  reaction-rale  profiles  of  sev¬ 
eral  important  branching-termination  reactions  for  diluted,  sto¬ 
ichiometric  methane-oxygen-nitrogen  flames  at  1  atm. 
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Fig.  14.  Numerically  calculated  reaction-rate  profiles  of  sev¬ 
eral  important  branching-termination  reactions  for  diluted,  sto¬ 
ichiometric  methane-oxygen-nitrogen  flames  at  4  atm. 


reaction  rate  profiles  for  the  three  dominant  chain 
reactions  for  the  84%-diluted  flames,  at  1  and  4 
atm,  respectively.  The  results  show  that  with  in¬ 
creasing  pressure  the  rates  of  the  termination  reac¬ 
tions  increase  faster  than  that  of  the  branching  re¬ 
action,  thus  further  substantiating  the  importance 
of  the  competition  between  the  branching  and  ter¬ 
mination  reactions  in  the  flame  behavior. 

CONCLUDING  REMARKS 

In  the  present  investigation  we  have  experi¬ 
mentally  determined  the  flame  speeds  of  nearly 
adiabatic,  stretch-free  methane-oxygen-nitrogen 
flames  that  are  subjected  to  variations  in  the  flame 
temperature  and  system  pressure.  The  viability  of 
these  data  is  further  substantiated  by  their  close 
agreement  with  the  numerically  calculated  values 
obtained  by  using  a  kinetic  scheme  previously  val¬ 
idated  for  the  determination  of  flame  speeds.  The 
availability  of  these  data  of  high  degree  of  accu¬ 
racy  not  only  is  required  for  the  present  study  but 
is  also  useful  for  future  modeling  efforts. 

Thus  by  analyzing  the  variations  of  the  overall 
reaction  order  and  activation  energy  with  pressure 
and  flame  temperature,  and  by  using  numerical 
sensitivity  analysis,  the  importance  of  termination 
reactions  in  the  overall  reaction  scheme  is  demon¬ 
strated  and  quantified.  The  influence  is  substan¬ 
tial  for  practically  all  of  the  situations  studied,  in¬ 
cluding  the  standard,  stoichiometric,  atmospheric 
pressure  case,  and  is  especially  significant  for  the 
weakly  burning,  highly  diluted,  elevated-pressure 


flames  that  are  relevant  to  flame  extinction  and 
stabilization  under  critical  conditions.  It  is  empha¬ 
sized  that  caution  is  needed  in  modeling  efforts 
that  use  a  one-step  overall  reaction  with  constant 
reaction  order  and  activation  energy,  because  the 
results  can  be  quantitatively  modified  and  some¬ 
times  even  qualitatively  affected  when  the  reaction 
order  changes  sign. 
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APPENDIX;  INFLUENCE  OF  BUOYANCY 
ON  FLAME  SPEED  DETERMINATION 

To  assess  the  influence  of  buoyancy  on  the  ac¬ 
curacy  of  the  counterflow  methodology  in  flame 
speed  determination,  especially  for  weak  llames, 
we  first  note  that  buoyancy  itself  does  not  di¬ 
rectly  affect  energy  and  species  transport  [19]. 
Indirectly,  it  modifies  the  flow  field  through  mo¬ 
mentum  transpon,  and  consequently  the  flame  re¬ 
sponse. 

Experimentally,  buoyancy  could  displace  the 
twin-flame  ensemble  upward  relative  to  the  syn; 
metry  plane.  However,  since  the  LDV  measure¬ 
ments  are  local  values,  the  modification  of  the 
flow  field  is  automatically  accounted  for  in  the 
measured  values.  To  substantiate  this  statement. 
Fig.  A1  shows  the  centerline  velocity  profile  of  a 
1-atm,  weakly  burning  (0  =  0.60),  CH^-air  flame 
ensemble.  It  is  seen  that  the  velocity  profiles  for 
the  upper  and  lower  flames  are  quite  symmetri¬ 
cal.  In  Fig.  A2  we  plot  the  reference,  stretched 
flame  speed  Sy  versus  the  stretch  rate  K  to  yield 
S°  zt  K  =0  through  extrapolation.  It  is  clear 
that  the  difference  in  the  determined  is  within 
the  range  of  experimental  uncertainty  stated  pre¬ 
viously. 

We  emphasize  that  the  key  ingredients  that  make 
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Fig.  A 1 .  Experimentally  determined  velocity  field  between  the 
two  nozzles  for  a  CFU-air  twm-flame  ensemble  with  p  -  1 
atm  and  4>  =  0.6. 

the  present  methodology  a  useful  one  are  (1)  the 
measurements  are  local  values,  and  (2)  extrap¬ 
olation  to  zero  stretch  eliminates  all  the  stretch 
effects,  whether  they  are  intentionally  imposed, 
or  inadvertently  present.  As  a  corollary,  we  cau¬ 
tion  against  the  use  of  the  bulk  stretch  rate,  ob¬ 
tained  f""  .formation  on  the  nozzle  separation 
dist  nec  aii'J  the  exit  velocity,  to  quantify  the  effect 
of  stretch  on  flames. 


K,  s  ‘ 


Fig  A2.  Reference,  stretched  flame  speed  as  a  function  of 
local  stretch  rate  for  the  upper  and  lower  flames 
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